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In this paper, we study the quantum properties of the three-mode squeezed operator. This 
operator is constructed from the optical parametric oscillator based on the three concurrent 



nonlinearities. We give a complete treatment for this operator including the symmetric 



and asymmetric nonlinearities cases. The action of the operator on the number and coherent 

Oh 

states are studied in the framework of squeezing, second-order correlation function, Cauchy- 
Schwartz inequality and single-mode quasiprobability function. The nonclassical effects are 
remarkable in all these quantities. We show that the nonclassical effects generated by the 

a -1 

asymmetric case-for certain values of the system parameters-are greater than those of the 
symmetric one. This reflects the important role for the asymmetry in the system. Moreover, 

> 

\Q ' the system can generate different types of the Schrodinger-cat states. 

- 

PACS numbers: 42.50.Dv,42.50.-p 



G\ '. I. INTRODUCTION 

O 



Squeezed light fulfils the uncertainty relation and has less noise than the coherent light in 
the one of the field quadratures. With the development of quantum information, squeezed states 
have become very important tool in providing efficient techniques for the encoding and decoding 
processes in the quantum cryptography For instance, the two-mode squeezed states are of great 
interest in the framework of continuous- variable protocol (CVP) [3]. In the CVP the quantum key 
distribution goes as follows. The two-mode squeezed source-such as parametric down conversion- 
emits two fields: one is distributed to Alice (A) and the other to Bob (B). Alice and Bob randomly 
choose to measure one of two conjugate field quadrature amplitudes. The correlation between the 
results of the same quadrature measurements on Alice's and Bob's side increases by increasing the 
values of the squeezing parameter. Through a public classical channel the users communicate their 
choices for the measurements. They keep only the results when both of them measure the same 
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quadrature and hence the key is generated. 

The generation of the multiparities squeezed entangled states is an essential issue in the multi- 



party communication including a quantum teleportation network [3|] , telecloning 



a, 



4j , and controlled 



dense coding [5]. The JV-mode CV entangled states have been generated by combining A' single- 
mode squeezed light in appropriately coupled beam splitters 3|. The three-mode CV entangled 
states have taken much interest in the literatures, e.g., 0, |f], 2, ls[ Q]- For instance, it has been 
theoretically shown that tripartite entanglement with different wavelengths can be generated by 
cascaded nonlinear interaction in an optical parametric oscillator cavity with parametric down 
conversion and sum- frequency generation [a]. Also, the three-mode CV states have been generated 



by three concurrent x^ nonlinearities [?|. This has been experimentally verified by the observa- 
tion of the triply coincident nonlinearities in periodically poled KT1OPO4 [8|]. More details about 
this issue will be given in the next section. Furthermore, the comparison between the tripartite 
entanglement in the three concurrent nonlinearities and in the three independent squeezed states 
mixed on the beam splitters , in the framework of the van Loock-Furusawa inequalities [l(J , has 
been performed in 7|- It is worth mentioning that the generation of macroscopic and spatially 
separated three-mode entangled light for triply coupled x Kerr coupler inside a pumped optical 
cavity has been discussed in 9]. It has been shown that the bright three-mode squeezing and full 
inseparable entanglement can be established inside and outside the cavity. 



Since the early days of the quantum optics squeezing is connected with what is so cal 



operator. There have been different forms of this operator in the literatures, e.g. 



15, n, H- 
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ed squeezed 
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For example, degenerate and non-degenerate parametric amplifiers are sources of the 
single- mode [11] and the two-mode [12] squeezing, respectively. The quantum properties of the 
three- mode squeezed operator (TMS), which is constructed from two parametric amplifiers and 
one frequency converter, have been demonstrated in [ijj]. This operator can be represented by the 



SU (1,1) Lie algebra generators 



]. Additionally, it can be generated-under certain condition- 
from bulk nonlinear crystal in which three dynamical modes are injected by three beams. Another 
possibility for the realization is the nonlinear directional coupler which is composed of two optical 
waveguides fabricated from some nonlinear material described by the quadratic susceptibility x ■ 
Finally, the mathematical treatments for particular type of the n-mode squeezed operator against 
vacuum states are given in 

In this paper, we treat the three concurrent parametric amplifiers given in [3, Q] as three- 
mode squeezed operator. We quantitatively investigate the nonclassical effects associated with this 
operator when acting on the three-mode coherent and number states. For these states we inves- 
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tigate squeezing, second-order correlation function, Cauchy-Schwartz inequality and single-mode 
quasiprobability functions. This investigation includes the symmetric (equal nonlinearities) and 
asymmetric (non-equal nonlinearities) cases. In the previous studies the entanglement of the sym- 



metric case only has been discussed 



The investigation in the current paper is motivated by 



the importance of the three concurrent parametric amplifiers in the quantum information research 



flfl 



] . Additionally, quantifying the nonclassical effects in the quantum systems is of fundamen- 
tal interest in its own right. We prepare the paper in the following order. In section 2 we construct 
the operator and write down its Bogoliubov transformations. In sections 3 and 4 we study the 
quadrature squeezing, the second-order correlation function as well as the Cauchy-Schwartz in- 
equality, respectively. In section 5 we investigate the single-mode quasiprobability functions. The 
main results are summarized in section 6. 

II. OPERATOR FORMALISM 

In this section we present the operator formalism for the optical parametric oscillator based on 
the three concurrent x^ nonlinearities. We follow the technique given in S] to construct the 
Hamiltonian of the system. In this regard, we consider three modes injected into a nonlinear crystal, 
whose susceptibility is x^ 2 \ to form three output beams at frequencies Uq, u)\, u)2- The interactions 
are selected to couple distinct polarizations. Assuming that x is the axis of the propagation within 
the crystal. The mode b\ is pumped at frequency and polarization (ujq + u>i,y) to produce the 
modes ai(uo,z) and 62(^1, y). The mode 62 is pumped at (ui\ + u>2,y) to produce the modes 0,2 
and 0,3 (u>2, z). Eventually, the mode 63 is pumped at (2wi, z) to produce the modes a\ and 0,2- The 
scheme for this interaction can be found in [3, Sj- The interaction Hamiltonian for this concurrent 



triple nonlinearity takes the form 



Hint = ih(xiha\al + X2&244 + X3^44) + n - c -> ( x ) 

where XjiJ = 1, 2, 3, represent the effective nonlinearities and h.c. stands for the hermitian conju- 
gate. The unitary operator associated with (P) is: 



U(t) = exp -it^* . (2) 
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In the undepleted pump approximation we set rj = Xj{bj(0)) as real parameters. Now we obtain 
the requested squeezed operator as: 

S(r) = exp[ri(aia 2 - a\a\) + r 2 (aia 3 - ajaj) + r 3 (a 2 a 3 - ajj&lj)], (3) 

where (r) = (r\,r2,r 3 ). Throughout this paper, the symmetric case means r\ = r<i = r 3 = r, 
otherwise we have an asymmetric case. It is evident that three disentangled state can be entangled 
under the action of this operator. This operator provides the following Bogoliubov transformations: 

St(r)o^(r) = f^ai + f®a\ + g?a 2 + g^a\ + hfa 3 + h^a\, j = 1, 2, 3 (4) 

where fjt,9ji,hjt,j > = 1,2 are functions in term of the parameters ri,r2,r 3 . The formulae of 
these functions for the asymmetric case are rather lengthy. Nevertheless, we write down only here 
the explicit forms for the symmetric case as Q]: 
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Relations (|3|) and © will be frequently used in the paper. For the symmetric case, the entanglement 
has been already studied in terms of the van Loock-Furusawa measure [7] . It has been shown that 
the larger the value of r, the greater the quantity of the entanglement in the tripartite. Moreover, 
the tripartite CV entangled state created tends towards GHZ state in the limit of infinite squeezing, 



3. 



In this paper we give an investigation for 



but is analogous to a W state for finite squeezing 
the entanglement of the asymmetric case from different point of view. This is based on the fact 
that the entanglement between different components in the system is a direct consequence of the 
occurrence of the nonclassical effects in their compound quantities and vice versa. We show that the 
asymmetric case can provide amounts of the nonclassical effects and / or entanglement greater than 
those of the symmetric case. Thus the asymmetry in the triply concurrent parametric amplifiers 
is important. 

The investigation of the operator ([3]) will be given through the three-mode squeezed coherent 
and number states having the forms: 

\4>n) = S(r)\ni,n 2 ,n 3 ), \ijj c ) = S(r)\a 1 ,a 2 ,a 3 ). (6) 
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Three-mode squeezed vacuum states can be obtained by simply setting rij = or ctj = in the 
above expressions. In the following sections we study the quantum properties for the states in 
greater details. 



Squeezing is an important phenomenon in the quantum theory, which can reflect the correlation 
in the compound systems very well. Precisely, squeezing can occur in combination of the quantum 
mechanical systems even if the single systems are not themselves squeezed. In this regard the 
nonclassicality of the system is a direct consequence of the entanglement. Squeezed light can 
be measured by the homodyne detector, in which the signal is superimposed on a strong coherent 
beam of the local oscillator. Additionally, squeezing has many applications in various areas, e.g., in 
quantum optics, optics communication, quantum information theory, etc [19]. Thus investigating 
squeezing for the quantum mechanical systems is an essential subject in the quantum theory. In 
this section we demonstrate different types of squeezing for the three-mode squeezed vacuum states 
([!]). To do so we define two quadratures X and Y, which denote the real (electric) and imaginary 
(magnetic) parts, respectively, of the radiation field as: 



where c\,C2 are c-numbers take the values or 1 to yield single-mode, two-mode and three-mode 
squeezing. These two operators, X and Y, satisfy the following commutation relation: 



where C = (1 + cf + c|)/2. It is said that the system is able to generate squeezing in the x- or 
y-quadrature if 



III. QUADRATURE SQUEEZING 




(7) 





(8) 
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C 



< 



X 



c 



or 



0) 




where 




Maximum squeezing occurs when S x 



= — 1 or 



Sy = "I. 
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For the symmetric case, one can easily deduce the following expressions: 
S x = 3(i+^|y {(l + c\ + C 2)[2exp(2r) + exp(-4r) - 3] + 2(ci + c 2 + cic 2 )[exp(-4r) - exp(2r)]}, 

Sy = 3(1+C i +C 2 ) {(l + c\ + C 2)[2exp(-2r) +exp(4r) - 3] + 2( Cl + c 2 + cic 2 )[exp(4r) - exp(-2r)]}. 



(10) 



For the single-mode case, c\ = c 2 = 0, the expressions (fTUj) reduce to: 

(11) 



5 X = i[2exp(2r) + exp(-4r) - 3], 



S y = |[2exp(-2r) + exp(4r) - 3]. 

It is evident that the system cannot generate single-mode squeezing. This fact is valid for the 
asymmetric case, too. For the two-mode case, c\ = l,c 2 = 0, i.e. first-second mode squeezing, we 
obtain 

S x = i[exp(2r) + 2exp(-4r) - 3], 
S y = ±[exp(-2r) + 2exp(4r) - 3]. 

Squeezing can be generated in the x-component only with a maximum value at r = ln(2)/3. Also 
the maximum squeezing, i.e. S x = — 1, cannot be established in this case. This is in contrast 
with the two-mode squeezed operator [12j for which S x = — 1 for large r. Roughly speaking, the 
quantum correlation in this system decreases the squeezing, which can be involved in the one of 
the bipartites. Finally, for three-mode case, c\ = c 2 = 1, we have 

S x = exp(-4r) - 1, S y = exp(4r) - 1. (13) 

Squeezing can be generated in the x-component only for r > 0. Squeezing reaches its maximum 
value for large values of r. The origin of the occurrence squeezing in (|13p is in the strong correlation 
among the components of the system. Moreover, the amount of the produced squeezing is two (four) 



12( (single-mode [111 ]) squeezed operator for certain values 



times greater than that of the two-mode 
of r. 

In Figs. [D[|a) and (b) we plot the squeezing parameter S x against r\ for two- and three-mode 
squeezing, respectively. We found that squeezing is not remarkable in S y . The solid curve is plotted 
for the symmetric case. The two-mode squeezing is given for the first-second mode system. We start 
the discussion with the two-mode case (Fig. [Ha)). From the solid curve, squeezing is gradually 
generated as n increases providing its maximum value S x = —0.206 at r = ln(2)/3 = 0.231, then 
it reduces smoothly and eventually vanishes, i.e. S x > 0, at r > ln(l + \/3)/2 = 0.5025. For the 
asymmetric case, squeezing increases gradually to be maximum S x = — 1 over a certain range of r±, 
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FIG. 1: Two- mode (a) and three- mode (b) squeezing against r\ for three- mode squeezed vacuum states. 
Solid, dashed and dotted curves are given for (Y 2 ,r 3 ) = (r"i,7"i), (0.1,0.2) and (0.4,0.6), respectively. 



then rapidly decreases and vanishes (see the dashed and dotted curves). This can be understood 
as follows. When the values of T2 and r% are relatively small, the main contribution in the system 
is related to the first parametric amplifier. Thus the system behaves as the conventional two-mode 
squeezed operator for a certain range of r±. This remark is noticeable when we compare the dotted 
curve to the dashed one. Generally, when the values of ri and increase, the degradation of the 
squeezing increases, too. Furthermore, in the range of r\ for which S x = — 1, the entanglement 
in the bipartite (1,2) is maximum, however, this is not the case for the other bipartities. This is 
connected with the fact: the quantum entanglement cannot be equally distributed among many 
different objects in the system. Comparison among different curves in Fig. 1(a) shows that the 
asymmetric case can provide amounts of squeezing much greater than those of the symmetric case. 
Now, we draw the attention to the three-mode squeezing, which is displayed in Fig. QJb). For 
the symmetric case, S x exhibits squeezing for r > 0, which monotonically increases providing 
maximum value for large r%, as we discussed above. This is in a good agreement with the fact 
that the symmetric case exhibits genuine tripartite entanglement for large values of r\ 21(. For 



the asymmetric case, the curves show initially squeezing, which reaches its maximum by increasing 
r± , then it gradually decreases and vanishes for large values of r\ . The greater the values of ri , r% 
the higher the values of the maximum squeezing in S x and the shorter the range of n over which 
squeezing occurs (compare dotted and dashed curves in Fig. 1(b)). This situation is the inverse 
of that of the two- mode case (compare Fig. [Ha) to (b)). Trivial remark, for small values of t\ the 
amounts of squeezing produced by the symmetric case are smaller than those by the asymmetric 
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one. We can conclude that for the asymmetric case the entanglement in the tripartite may be 
destroyed for large 7*1, where S x > 0. Of course this is sensitive to the values of r2,r^. Conversely, 
the amounts of the entanglement between different bipartites in the system for the asymmetric 
case can be much greater than those of the symmetric one for particular choice of the parameters. 
The final remark, the amounts of squeezing produced by the operator ([3]) are greater than those 
generated by the TMS Q. 



IV. SECOND-ORDER CORRELATION FUNCTION AND CAUCHY-SCHWARTZ 

INEQUALITY 

In this section we study the second-order correlation function and the Cauchy- Schwartz inequal- 
ity for the states ©. These two quantities are useful for getting information on the correlations 
between different components in the system. In contrast to the quadrature squeezing, these quan- 
tities are not phase dependent and are therefore related to the particle nature of the field. We start 
with the single-mode second-order correlation function, which for the jth mode is defined as: 

S f'(0) = ^t-1, (14) 



(2) (2) 

where g\ (0) = for Poissonian statistics (standard case), g)- (0) < for sub-Poissonian statistics 

(2) 

(nonclassical effects) and g- (0) > for super-Poissonian statistics (classical effects). The second- 
order correlation function can be measured by a set of two detectors 22], e.g. the standard Hanbury 



BrownTwiss coincidence arrangement. Furthermore, the sub-Poissonian light has been realized in 



the resonance fluorescence from a two- level atom driven by a resonant laser field [23|. We have 
found that three-mode squeezed coherent state cannot yield sub-Poissonian statistics. Thus we 
strict the study here to g[ (0) of the three-mode squeezed number states. Form (jlj) and pj), one 
can obtain the following moments: 
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(a) 




-0.5 



(b) 



FIG. 2: Second-order correlation of the first mode against r\. (a) (711,712,713) = (1,1,1), (r2,rs) = {r\,ri) 
solid curve, (0.1,0.2) dashed curve, and (0.4,0.6) dotted curve, (b) (r- 2 ,r 3 ) = (rx,rx), (n±, n 2 , n 3 ) = (1,0,0) 
solid curve and (0, 1, 0) dashed curve. 



(a[ai) = mff>' + (m + l)ff>* + n 2 g\^ + (n 2 + 1)^ + n 3 /7,^ + (n 3 + 1)^, 



f (l)2 



f (l)2 



,(1)2 



(1)2 



,(1)2 



,(1)2 



(a\ 2 al) = ni ( ni - l)/^ 4 + (ni + l)(m + 2)/ 2 W4 + (2m + l) 2 /} J / 2 l 



+n 2 (n 2 - 1) 5 [ 1)4 + (n 2 + l)(n 2 + 2)g 2 ^ + (2n 2 + l)' 2 g^g y 2 



(1)4 



, 2 (1)2 (1)2 



+n 3 (n 3 - l)h[ 1)4 + (n 3 + l)(ra 3 + 2)h^ + (2n 3 + l) 2 h^ )2 h { 2 



+(2m + l)/ 1 w / 2 w [2(2n a + l)g? } gr + (2n 3 + l)h\ L % 



(15) 



+(2n 3 + l)^ 1) 4 1) [2(2n 2 + 1)^^^^ + (2m + 1)/^ 



(i)„(i) 



r(l) f(l)l 



+4[n 1 / 1 {1)2 + (ni + l)/^Hhsi^ + in 2 + + n 3 h[ 1,z + (n 3 + 



,(1)2! r (1)2 



(1)2 



,(1)2 



,(1) 2 1 



+4[n 25 f )2 + (n 2 + l)^ 1)2 ][n 3 ^ 1)2 + (n 3 + 1)4 1)2 ]. 

By means of (j!4j) and (|15p . the quantity (0) is depicted in Figs. [2] for given values of the 

parameters. In Fig. [5]^a) we present the role of the squeezing parameters r* on the behavior of 

(2) 

the g{ (0). From the solid curve (, i.e. symmetric case), one can observe that the maximum sub- 
Piossonian statistics occur for relatively small values of rj. In this case, the system tends to the 
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Fock state |1), which is a pure nonclassical state. As the values of r\ increase, the nonclassicality 
monotonically decreases and completely vanishes around n — 0.3. Comparison among the curves 
in Fig. 2(a) shows when the values of (r2, r^) increase, the amounts of the sub-Poissonian statistics 
inherited in the first mode decrease. This is connected with the nature of the operator ([3|), in 
which the behavior of the single-mode undergoes an amplification process caused by the various 
down-conversions involved in the system. Furthermore, particular values of the asymmetry can 
enlarge the range of nonclassicality (compare the solid curve to the dashed one). Now, we draw 
the attention to the Fig. 2(b), which is given to the symmetric case . From this figure one realizes 
how can obtain sub-Poissonian statistics from a particular mode as an output from the operator 
(|3|). Precisely, this mode should be initially prepared in the nonclassical state. We can analytically 
prove this fact by substituting n\ = 0, rt2 = = n into (fT4"|) and (fT5|) , After minor algebra, we 
arrive at: 

{afa\) - (ajai) 2 = / 2 (1)4 + 2n(n - l)gf )A + 2(n + l)(n + 2)^ 1)4 

+2n(n + 1)^ 1)2 4 1)2 + [f?^ + 2(2n + l)^] 2 (16) 
+4/( 1)2 [n^ 1)2 + (n + l)^ 1)2 ]>0. 

The final remark, the comparison between the solid curves in Figs. 2(a) and (b) shows that the 
nonclassical range of r\ in (b) is greater than that in (a). In other words, to enhance the sub- 
Poissonian statistics in a certain mode, the other modes have to be prepared in states close to the 
classical ones. 




(a) (b) (c) 

FIG. 3: The parameter Vjk against r\ for the coherent state with oi\ — 0,2 = «3 = 1, where (j, k) = (1, 2) 
(a), (1,3) (b) and (2, 3) (c). Additionally, (r 2 , r 3 ) = (n, ri) solid curve, (0.1, 0.2) dashed curve and (0.4, 0.6) 
dotted curve. 



The violation of the classical inequalities has verified the quantum theory. Among these in- 
equalities is the Cauchy-Schwarz inequality 24|, which its violation provides information on the 
intermodal correlations in the system. The first observation of this violation was obtained by 
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Clauser, who used an atomic two-photon cascade system [25(1 . More recently, strong violations us- 
ing four- wave mixing have been adopted in H]- I* 1 addition, a frequency analysis has been used 



to infer the violation of this inequality over a limited frequency regime 
inequality is Vjk < 0, where Vjk has the form: 



28[. The Cauchy-Schwarz 



(17) 



Occurrence of the negative values in Vjk means that the intermodal correlation is larger than the 
correlation among the photons in the same mode. This indicates a strong deviation from the 
classical Cauchy-Schwarz inequality. This is related to the quantum mechanical features, which 
include pseudodistributions instead of the true ones. In this respect, the Glauber-Sudarshan P 
function possesses strong quantum properties [3]. 




FIG. 4: The parameter Vjk against n for the Fock-state case with n\=n2 = TI3 = 1, where (j, k) = (1, 2). 
Additionally, (r2,r 3 ) = (Vi,ri) solid curve, (0.1,0.2) dashed curve and (0.4,0.6) dotted curve. 



We have found that the Cauchy-Schwartz inequality can be violated for both coherent- and 
number-state cases. The expressions of the different quantities in (|17p are too lengthy but straight- 
forward and hence we don't present them here. We start with the three-mode squeezed coherent 
states. Information about them is shown in Figs. 3(a), (b) and (c), for given values of the system 
parameters. The negative values are remarkable in most of the curves, reflecting the deviation 
from the classical inequality. For the symmetric case the nonclassical correlation is remarkable 
for n > 0.4, which increases gradually till r\ ~ 0.6 providing its maximum value, then smoothly 
reduces and eventually goes to zero, i.e. Vjk — 0, for large r\. For the asymmetric case, the 
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deviation from the classical inequality is obvious, which may be smaller or greater than those in 
the symmetric one based on the competition among different nonlinearities in the system rj. In 
other words, e.g., the correlation between modes 1 and 2 is much stronger than that between the 
others only when r\ > r2,r%. This can be understood from the structure of the operator ([3]). 
Comparing this behavior to that of the two-mode squeezing given in the preceding section one can 
conclude that a large amount of squeezing does not imply large violation of the inequality [30] . As 
we mentioned before: the entanglement is a direct consequence of the occurrence of the nonclas- 
sical effects. As a result of this, the behavior of the two-mode squeezing and Vjk may provide a 
type of contradiction. Precisely, the bipartite can be entangled (non-entangled) with respect to, 
say, two-mode squeezing (Vjk)- This supports the fact that these two quantities provide only a 
sufficient condition for entanglement. Considering both of them we may obtain conditions closer 
to the necessary and sufficient condition. A study about this controversial issue has been already 
discussed for the entanglement in a parametric converter [20], where different entanglement criteria 
leaded to different results. 

In Fig. 4 we plot the parameter Vjk for the Fock-state case. From this figure the deviation 
from the classical Vjk is quite remarkable. For the symmetric case, maximum deviation occurs in 
V\ t 2 for r = 0, monotonically decreases as r evolving and vanishes at r ~ 1. Comparison among 
different curves in this figure shows that the asymmetry can enlarge the range of t\ over which 
the deviation of the inequality occurs. Similar behavior has been observed for V2 3 and Vjj. 3 (we 
have checked this fact). In conclusion, the violation of the classical inequalities provides an explicit 
evidence of the quantum nature of intermodal correlation between modes. This is not surprising, 
as the entanglement is a pure quantum mechanical phenomenon that requires a certain degree of 
nonclassicality either in the initial state or in the process that governs the system. 



V. QUASIPROB ABILITY DISTRIBUTION FUNCTION 



Quasiprobability distribution functions, namely, Husimi function (Q), Wigner function (W), and 
Glauber P functions, are very important since they can give a global description of the nonclassical 
effects in the quantum systems. These functions can be measured by various means, e.g. photon 
counting experiments [311 ]. using simple experiments similar to that used in the cavity (QED) and 



ion traps 



32 



33j, and homodyne tomography For the system under consideration, we focus 



the attention here on the single-mode case, say, the first mode for the three-mode squeezed number 
states 0. We start with the s— parameterized characteristic function C(C,s), which is defined as: 
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C(C, s) = Tr[pex V ((a\ - C&i + ^\C\% (18) 

where p is the density matrix of the system under consideration and s is a parameter taking the 
values 0, 1, —1 corresponding to symmetrically, normally and antinormally ordered characteristic 
functions, respectively. For three-mode squeezed number states and with the help of the relations 
([3]) one can easily obtain: 



C(C, s) = exp[-iM 2 - in 2 - i|t, 3 | 2 + 1 |CI 2 ] 
xL ni (M 2 )L n2 (M 2 )L n3 (M 2 ), 



(19) 



where 



v 1 = Cfl 1) -Cti\ - 2 = rf ) -C^ 1) , v 3 = C h?-Ch^ (20) 

and L^(.) is the associated Laguerre polynomial having the form: 

L ^-^ oil+mk -i)iiv (21) 

The s-parameterized quasiprobability distribution functions are defined as 

W(z, s) = it- 2 J d 2 (C((, s) exp« - C**)> ( 22 ) 

where z = x + iy and s = 0, 1, —1 are corresponding to W, P, Q functions, respectively. On substi- 
tuting (|19p into (|22p and applying the method of the differentiation under the sign of integration 
we can obtain the following expression: 

{ni,n 2 ,n 3 } {j',j+k} I n 

w(z,s) = ± E E f ' 



j',j,k=0 h,h=0 \ j' 




j'.j'lk 



v / f (l)2 , -(l)2uw (1)2 (l)2w 2 /,(l),(l)^_j / (1) (l)w +fc _, gh+h, Qi+S+k-h-h | ( 23 ) 



where 



A 1 = / 1 (1)2 + /p +5 i 1)2 + ^ 1)2 + / i f )2 + 4 1)2 , 



A 2 = / 1 (1) / 2 (1) +^ 1) 4 1) + ^4 1) , (24) 



B = i(Ax - s) - h, K = B 2 -{A 2 + b 2 ) 2 . 
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The correlation between the modes in the system can be realized in W(z, s) as cross terms, e.g. in 
A2. This can give a qualitative information about the entanglement in the system. For the three- 
mode squeezed vacuum states (, i.e., n\ = 77,2 = 77.3 = 0) the W function (|23l) can be expressed 



as: 



1 x 2 y 2 



where 



W( X ,y,s) = - 7 —eM-^-^l (25) 



^ + = 2((AA > 1 ) 2 )-§, 

= |[(/ 1 ( n/f) 2 + (^ + ^) 2 +(4 1) +4 1) ) 2 ]-i, 

= 2((AY 1 ) 2 ) - f , 



(26) 



= |[(/ 1 (1) -/i 1) ) 2 + (^ 1) -^ 1) ) 2 + (4 1) -4 1) ) 2 ]-f- 

From (|25p and (|26p it is evident that the quasidistributions are Gaussians, narrowed in the y 
direction and expanded in the x direction. Nevertheless, this does not mean squeezing is available 
in this mode. Actually, this behavior represents the thermal squeezed light, which, in this case, 
is a super-classical light. Precisely, with s = the W function exhibits stretched contour, whose 
area is broader than that of the coherent light. In this regard the phase distribution and the 
photon-number distribution associated with the single-mode case exhibits a single-peak structure 
for all values of rj. This peak is broader than that of the coherent state, which has the same mean- 



jhoton number. Actual 
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16, 



y, this is a quite common property for the multimode squeezed operators 



17(. Thus the single-mode vacuum or coherent states, as outputs from 
the three concurrent amplifiers described by ([3]) are not nonclassical states. This agrees with the 
information given in the sections 3 and 4. 

Now we consider two cases: (711,77.2,713) = (0,0,77,3) and (77-1,77,2,773) = (77,1, 0, 0). For the first 
case, we investigate the influence of the nonclassicality in the third mode on the behavior of the 



first mode. To do so we substitute n\ = n<i = in (|23p and after minor algebra we arrive at: 



(27) 



E 

771=0 
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where 

V± = (h?±h^f-^. (28) 

One can easily check when r\ = T2 = r% = the W function (|27|) reduces to that of the vacuum 
state. The form (]27p includes Laguerre polynomial, which is well known in the literature by 
providing nonclassical effects in the phase space. This indicates that the nonclasssical effects can 
be transferred from one mode to another under the action of the operator ([3]). Of course the 
amount of the transferred data depends on the values of the squeezing parameters rj . 




(a) 



o 
o 




1.5 2 2.5 3 3.5 



(b) 



FIG. 5: The Wigner function of the first mode (a) and the evolution of the phase space origin of the Wigner 
function (b). In (a) we use {r\, r2, r^, n\, n^, n^) = (1.1,1.1,1.1,0,0,1). In (b) we use (r±, r^, n\, 122, n^) = 
(r 3 , r 3 , 0, 0, 1) solid curve, (0.6, 0.8, 0, 0, 1) dashed curve, and (0.6, 0.8, 1, 0, 0) dotted curve. 



The second case (ni,0,0) has the same expression (j27|) with the following transformations: 



(1) , ,(lh2 



n 3 -> ni, r]± = (f± ± f. 



(29) 



Now we prove that this W function tends to that of the number state when r\ = r2 = r% = 0. In 
this case, the transformations (1291) tend to: 



1+a 1-8 

V± = — ^— > #± = — ^— • 



(30) 
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(a) 



(b) 





(c) 




(d) 



FIG. 6: The Wigner function of the first mode when (r±, Ti, r^, n\, na, TI3) 
(0.6,0.8,0.9,0,0,1) (b), (0.6,0.8,2,0,0,1) (c) and (0.4,0.8,2,0,0,2) (d). 



(0.6,0.8,0.9,1,0,0) (a), 



Substituting these variables in the expression (|27p (with 713 — > ni) we obtain: 

2(-i)"i (l+s)"^ exp ^ 2{x*+y*y 2i 

m=0 



' , PTn [ 2(x 2 +y 2 ) 1 t ~f r 4a: 2 it ~f r iy" j 



_ 2(-l)*i r 2(x 2 +y 2 ) lT r 4(x 2 



(31) 



7T (1- S ) 



— exp[ — \h ni [ 1 _ s2 j, 



In (|3ip the transition from the first line to the second one has been done using the identity: 

m 

2>£(x)L£_ n (j,) = L^ 2 +V + y)- (32) 



n=0 



The expression (I3ip is the s-quasprobability distribution for the number state, e.g. 



35J. We 
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conclude this part by writing down the form of the W function at the phase space origin, which 
is a sensitive point for the occurrence of the nonclassical effects. Moreover, it can simply give 
visualization about the behavior of the system. Additionally, this point can be measured by the 
photon counting method [31]. From (|27p we have: 



It is obvious that the Winger function exhibits negative values at the phase space origin only when 
713 is an odd number. 

In Figs. 5 and 6 we plot the W functions for the given values of the system parameters. We start 
the discussion with the symmetric case. In Fig. 5(a) we use (ni, 712,71-3) = (0,0,1) meaning that 
the mode under consideration is in the vacuum state. Thus for rj = the W function exhibits the 
single-peak-Gaussian structure with a center at the phase space origin. When rj 7^ this behavior 
is completely changed, where one can observe a lot of the nonclassical features, e.g. negative values, 
multipeak structure and stretching contour (see Fig. 5(a)). This indicates that the nonclassical 
effects can be transferred from one mode to the other under the action of the operator ([3]). In 
Fig. 5(b) we plot the "evolution" of the W function given by f|33|) against the parameter r (7*3) 
for the symmetric (asymmetric) case. The aim of this figure is to estimate the exact value of 
the nonlinearity r (7-3) for which the nonclassical effects maximally occur and/or transfer from 
certain mode to the other. For the symmetric case, this occurs at r = 1.2, while the nonclassicality 
is completely washed out at r = 3. Now, we draw the attention to the asymmetric case which is 
plotted in Figs. 6. Fig. 6(a) gives information on the case (711,712,713) = (1,0,0). The W function 
of the Fock state |1) is well known in the literatures by having inverted peak in phase space with 
maximum negative values. This is related to that this state provides maximum sub-Poissonian 
statistics. Under the action of the operator ([3]) these negative values are reduced and the two- 
peak structure is started to be constructed. This indicates that the system is able to generate 
particular types of the Schrddinger-cat states by controlling the system parameters. This is really 
obvious in Figs. 6(b)-(d), which are given to the cases 711=712 = 0,713 = 1,2. For instance, from 
Fig. 6(b) the W function provides a two-peak structure. Nevertheless, by increasing the values 
of the parameter 7*3, the W function exhibits the two Gaussian peaks and inverted negative peak 
in-between indicating the occurrence of the interference in phase space (see Fig. 6(c)). This shape 
is similar to that of the odd-coherent state. Additionally, the Fig. 6(d), in which 713 = 2, provides 
the well-known shape of the W function of the even coherent state. Generally, the even and the 
odd Schrodinger-cat states have nearly identical classical components (, i.e. the positive peaks) 




(33) 
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and only differ in the sign of their quantum interferences. These are interesting results, which 
show that by controlling the nonlinearity of the system and preparing a certain mode in the Fock 
state |1) or |2) one can generate cat states. It is worth mentioning that the Fock state \n) can be 



37 L Similar results have 



36j |. Furthermore, quite 



prepared with very high efficiency according to the recent experiments 
been obtained from the codirectional three-mode Kerr nonlinear coupler 
recently the construction of the cat state trapped in the cavity in which several photons survive 
long enough to be repeatedly measured is given in 



. In this technique, the atoms crossing the 



cavity one by one are used to obtain information about the field. We proceed, we have noted that 
the W function of the cases ri2 = n.3 = 0,ni = 1,2 can provide quite similar behaviors as those 
of Figs. 6(c) and (d) for the same values of rj. Now, we draw the attention to the dashed and 
dotted curves in the Fig. 5(b). These curves provide information on the evolution of the VF(0,0) 
against r% for the case of Fig. 6(a) and (c), respectively. From the dotted curve, i.e. the mode 
under consideration is in the Fock state, the W(0, 0) exhibits the maximum negativity at r% = 0, 
which monotonically decreases and completely vanishes at = 4. This shows for how long the 
nonclassicality inherited in the first mode survives based on the intensity of the third amplifier. 
On the other hand, form the dashed curve, i.e. the mode under consideration is in the vacuum 
state, the W(0, 0) exhibits negative values for r% > 1, increases rapidly to show maximum around 
r3 = 2, reduces gradually and vanishes for r% > 5. This range of negativity is greater than that 
of the dotted curve. This is in a good agrement with the behavior of the second-order correlation 
function. Finally, comparison among the curves in Fig. 5(b) confirms the fact: for certain values 
of the system parameters, the asymmetric case can provide nonclassical effects greater than those 
of the symmetric one. 



VI. CONCLUSION 



In this paper we have studied the three-mode squeezed operator, which can be implemented 
from the triply coincident nonlinearities in periodically poled KTiOPO^. The action of this oper- 
ator on the three-mode coherent and number states is demonstrated. We have studied quadrature 
squeezing, second-order correlation function, Cauchy- Schwartz inequality and quasiprobability dis- 
tribution function. The obtained results can be summarized as follows. Generally, the single-mode 
vacuum or coherent states, as outputs from the three concurrent amplifiers are not nonclassical 
states. The system can exhibit two-mode and three-mode squeezing. The amount of the two-mode 
squeezing generated by the asymmetric case is much greater than that of the symmetric case. 
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Three-mode squeezed coherent (number) states cannot (can) exhibit sub-Poissonian statistics. To 
obtain maximum sub-Poissonian statistics from a particular mode, under the action of the op- 
erator ([3]), it must be prepared in the nonclassical state and the other modes in states close to 
the classical ones. We have found that the Cauchy- Schwartz inequality can be violated for both 
coherent states and number states. The origin in the violation is in the strong quantum correla- 
tion among different modes. For the Fock-state case, the asymmetry in the system enhances the 
range of nonlinearities for which Vjk is nonclassical compared to that of the symmetric one. In 
the framework of the quasiprobability distribution we have shown that the nonclassical effects can 
be transferred from one mode to another under the action of the operator ((3|). The amount of 
transferred nonclassicality is sensitive to the values of the squeezing parameters. Interestingly, the 
system can generate particular types of the Schrodinger-cat states for certain values of the system 
parameters. Generally, we have found that the nonclassical effects generated by the operator ([3]) 
are greater than those obtained from the operator TMS [14]. Finally, the asymmetry in the three 
concurrent nonlinearities process is important for obtaining significant nonclassical effects. 
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